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This  paper  reviews  the  development  of  latent  heat  thermal  energy  storage  systems  studied  detailing 
various  phase  change  materials  (PCMs)  investigated  over  the  last  three  decades,  the  heat  transfer  and 
enhancement  techniques  employed  in  PCMs  to  effectively  charge  and  discharge  latent  heat  energy  and 
the  formulation  of  the  phase  change  problem.  It  also  examines  the  geometry  and  configurations  of  PCM 
containers  and  a  series  of  numerical  and  experimental  tests  undertaken  to  assess  the  effects  of 
parameters  such  as  the  inlet  temperature  and  the  mass  flow  rate  of  the  heat  transfer  fluid  (HTF).  It  is 
concluded  that  most  of  the  phase  change  problems  have  been  carried  out  at  temperature  ranges  between 
0  °C  and  60  °C  suitable  for  domestic  heating  applications.  In  terms  of  problem  formulation,  the  common 
approach  has  been  the  use  of  enthalpy  formulation.  Heat  transfer  in  the  phase  change  problem  was 
previously  formulated  using  pure  conduction  approach  but  the  problem  has  moved  to  a  different  level  of 
complexity  with  added  convection  in  the  melt  being  accounted  for.  There  is  no  standard  method  (such  as 
British  Standards  or  EU  standards)  developed  to  test  for  PCMs,  making  it  difficult  for  comparison  to  be 
made  to  assess  the  suitability  of  PCMs  to  particular  applications.  A  unified  platform  such  as  British 
Standards,  EU  standards  needs  to  be  developed  to  ensure  same  or  similar  procedure  and  analysis 
(performance  curves)  to  allow  comparison  and  knowledge  gained  from  one  test  to  be  applied  to  another. 
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1.  Introduction 

The  quest  for  new  technologies  to  avert  the  growing  concern 
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Nomenclature 

A  area  (m2);  thermal  diffusivity  ratio,  A  =  aslotf 

a  equation  coefficient 

Bi  Biot  number,  Bi  =  hrQIX 

B  solid  subcooling  parameter,  B  =  XSA Ts/X  fATsf 

b  width  of  a  rectangular  PCM  container  (m) 

C  instantaneous  height-averaged  cooling  rate 

Cp  specific  heat  capacity  (kj/kg  K) 

d  diameter  (m) 

Fo  Fourier  number 

/  liquid  fraction  of  PCM  melt 

g  acceleration  due  to  gravity  (m/s) 

Gr  Grashof  number 

H  enthalpy  (J) 

h  heat  transfer  coefficient  (W/(m2  K)) 

He  instantaneous,  height-averaged  melting  rate 

HTF  heat  transfer  fluid 

k  thermal  conductivity  (W/m) 

l  length  (m) 

M  instantaneous,  height-averaged  melting  rate 

m  mass  (kg) 

rnmass  flow  rate  (kg/s) 

Pr  Prandtl  number 

PCM  phase  change  material 

Q.  energy  or  heat  flow  (J) 

q  heat  flux  (W/m2) 

R  thermal  resistance  (K/W) 

Ra  Rayleigh  number 

Re  Reynold’s  number 

r  radius  (m) 

Ste  Stefan  number  or  liquid  superheating,  Ste  =  CATfi 
XT, 

T  temperature  (°C  or  I<) 

AT  temperature  difference  (°C  or  I<) 

t  time  (s) 

At  time  step  (s) 

v  velocity  components  (m/s) 

V  volume  (m3) 

v  velocity  (m/s) 

Xf  local  volume  fraction  of  fibers 

x,  y,  z  x,  y,  z-component  in  a  Cartesian  coordinate  system 

Ax  space  increment  (m) 

w  dimensionless  geometric  parameter  defined  by  1  + 

w  =  lA/V 

Z  fin  height  (m) 

Greek  symbols 

8  fin  thickness  (m);  dimensionless  front  radius 

(interface  position) 

a  thermal  diffusivity  of  PCM  (m2/s) 

coefficient  of  volumetric  expansion  (K_1) 
p  density  (kg/m3) 

A  latent  heat  of  fusion  (kj/kg) 

6  angle  (°) 

8  dimensionless  front  radius,  8  =  r/rG 

O  dimensionless  temperature  difference 

v  dynamic  viscosity  (m2/s) 

pt  viscosity  (kg/(m  s)) 


8  emissivity,  emmittance 

r  dimensionless  time 


Subscripts  and  superscripts 


c 

convective 

Cd 

conductive 

c 

cylindrical 

cr 

critical 

e 

enclosure 

eff 

effective 

eq 

equivalent 

f 

fluid 

fiber 

fiber 

fin 

fin 

fric 

friction 

H 

heat  transfer  fluid 

i 

inner  (inlet) 

l 

liquid 

lag 

lag 

m 

melting 

max 

maximum 

0 

outer  (outlet) 

P 

PCM 

r 

rectangular  container 

ref 

reference 

s 

solid 

t 

tube 

concern  over  the  last  three  decades.  Central  to  the  problem  is  the 
need  to  store  excess  energy  that  would  otherwise  be  wasted  and 
also  to  bridge  the  gap  between  energy  generation  and  consump¬ 
tion.  Latent  heat  thermal  energy  storage  is  particularly  attractive 
technique  because  it  provides  a  high  energy  storage  density.  When 
compared  to  a  conventional  sensible  heat  energy  storage  systems, 
latent  heat  energy  storage  system  requires  a  smaller  weight  and 
volume  of  material  for  a  given  amount  of  energy.  In  addition  latent 
heat  storage  has  the  capacity  to  store  heat  of  fusion  at  a  constant  or 
near  constant  temperature  which  correspond  to  the  phase 
transition  temperature  of  the  phase  change  material  (PCM).  The 
study  of  phase  change  materials  was  pioneered  by  Telkes  and 
Raymond  [1  ]  in  the  1940s  but  did  not  receive  much  attention  until 
the  energy  crisis  of  late  1970s  and  early  1980s  where  it  was 
extensively  researched  for  use  in  different  applications  especially 
for  solar  heating  systems  [2-33  ].  Less  attention  was  however  given 
to  phase  change  storage  research  after  the  energy  crisis  until  it 
emerged  again  recently.  Since  the  late  70s  a  number  of  studies 
have  been  conducted  to  assess  the  overall  thermal  behaviours  of 
latent  heat  thermal  storage  systems  [34-59].  Studies  of  phase 
change  systems  have  investigated  design  fundamentals,  system 
and  process  optimization,  transient  behaviour,  and  field  perfor¬ 
mance.  The  research  and  development  has  been  broad  based  and 
productive,  concentrating  on  both  the  resolution  of  specific  phase 
change  materials  and  problems  and  the  study  of  the  characteristics 
of  new  materials.  The  major  disadvantage,  as  reported  by  many 
researchers  has  been  the  low  thermal  conductivities  possessed  by 
many  PCMs,  leading  to  low  charging  and  discharging  rates 
(especially  for  the  organic  based  materials).  The  development  of 
a  latent  heat  thermal  energy  storage  system  therefore  involves  the 
understanding  of  heat  transfers/exchanges  in  the  PCMs  when  they 
undergo  solid-to-liquid  phase  transition  in  the  required  operating 
temperature  range,  the  design  of  the  container  for  holding  the  PCM 
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and  formulation  of  the  phase  change  problem.  This  paper  reviews 
the  studies  conducted  on  phase  change  materials  for  different 
applications  and  reported  in  the  open  literature. 

2.  Phase  change  materials  (PCMs) 

Energy  storage  may  be  in  the  form  of  sensible  heat  in  a  liquid  or 
solid  medium,  as  heat  of  fusion  (latent  heat),  or  as  chemical  energy 
or  products  in  a  reversible  chemical  reaction.  The  classification  of 
energy  storage  and  the  materials  used  are  detailed  in  [51]. 
Chemical  energy  storage  has  not  as  yet  been  used  in  practical 
applications  and  both  technical  and  economical  questions  have  yet 
to  be  answered  for  some  of  the  possibilities  proposed.  To  date  most 
of  the  studies  conducted  on  storage  materials  have  concentrated 
on  sensible  and  latent  heat  storage  systems.  Studies  conducted  to 
compare  phase  change  and  sensible  heat  storages  have  shown  that 
a  significant  reduction  in  storage  volume  can  be  achieved  using 
PCM  compared  to  sensible  heat  storage.  Results  of  studies  from 
Morrison  and  Abdel-Khalik  [48]  and  Ghoneim  [60]  show  that  to 
store  the  same  amount  of  energy  from  a  unit  collector  area,  rock 
(sensible  heat  storage  material)  requires  more  than  seven  times 
the  storage  mass  of  Paraffin  116  Wax  (PI  16-Wax),  five  times  the 
storage  mass  of  medicinal  paraffin  and  more  than  eight  times  the 
storage  mass  of  Na2SO410H2O. 

2.2.  PCMs  investigated 

Several  authors  have  carried  out  investigation  into  a  wide  range 
of  PCMs,  subdividing  them  into  organic,  inorganic,  eutectics  PCMs 
[61].  The  main  criteria  that  govern  the  selection  of  phase  change 
heat  storage  materials  are  [51]: 

•  Possess  a  melting  point  in  the  desired  operating  temperature 
range  (temperature  range  of  application). 

•  Possess  high  latent  heat  of  fusion  per  unit  mass,  so  that  a  smaller 
amount  of  material  stores  a  given  amount  of  energy. 

•  High  specific  heat  to  provide  additional  significant  sensible  heat 
storage  effects. 

•  High  thermal  conductivity,  so  that  the  temperature  gradients  for 
charging  and  discharging  the  storage  material  are  small. 

•  Small  volume  changes  during  phase  transition,  so  that  a  simple 
container  and  heat  exchanger  geometry  can  be  used. 

•  Exhibit  little  or  no  subcooling  during  freezing. 

•  Possess  chemical  stability,  no  chemical  decomposition  and 
corrosion  resistance  to  construction  materials. 

•  Contain  non-poisonous,  non-flammable  and  non-explosive 
elements/compounds. 

•  Available  in  large  quantities  at  low  cost. 

Table  1  lists  the  current  companies  that  commercially  produce  over 
100  PCMs.  In  addition  to  these,  several  PCMs  have  been  proposed  or 
studied  by  different  researchers.  A  detailed  list  of  PCMs  studied  or 
proposed  for  study  can  be  found  in  [51,62-65].  Table  2  gives  the 


Table  1 

Commercial  PCM  manufacturers  in  the  world. 


Manufacturer 

PCM  temperature 
range 

Number  of 
PCMs  listed 

RUBITHERM  (www.rubitherm.de) 

-3  °C  to  100°C 

29 

Cristopia  (www.cristopia.com) 

-33  °C  to  27  °C 

12 

TEAP  (www.teappcm.com) 

-50  °C  to  78  °C 

22 

Doerken  (www.doerken.de) 

-22  °C  to  28  °C 

2 

Mitsubishi  Chemical  (www.mfc.co.jp) 

9.5  °C  to  118  °C 

6 

Climator  (www.climator.com) 

-18  °C  to  70  °C 

9 

EPS  Ltd  (epsltd.co.uk) 

—  1 14 °C  to  164°C 

61 

Merck 

thermophysical  data  for  a  range  of  typical  and  important  PCMs 
studied  or  proposed  for  study  by  previous  authors. 

2.2.  Classification  of  PCMs  to  melting  temperature  range  and 
application  area 

The  selection  of  an  appropriate  PCM  for  any  application 
requires  the  PCM  to  have  melting  temperature  within  the  practical 
range  of  application.  Several  application  areas  have  been  proposed 
for  PCMs  studied.  Table  2  lists  some  of  the  target  application  areas 
for  the  selection  of  PCMs  for  study,  as  provided  by  the  references  in 
the  open  literature.  From  Table  2,  it  can  be  seen  that  most  of  the 
research  on  phase  change  problems  have  been  carried  out  within 
the  temperature  range  0-65  °C  suitable  for  domestic  heating/ 
cooling  (Table  3). 

3.  PCM  containers 

Once  the  PCM  has  been  selected  based  primarily  on  the 
temperature  range  of  application,  the  next  most  important  factors 
to  consider  are: 

(i)  the  geometry  of  the  PCM  container  and 

(ii)  the  thermal  and  geometric  parameters  of  the  container 
required  for  a  given  amount  of  PCM. 

Each  of  these  factors  has  a  direct  influence  on  the  heat  transfer 
characteristics  in  the  PCM  and  ultimately  affects  the  melt  time  and 
the  performance  of  the  PCM  storage  unit. 

To  ensure  long-term  thermal  performance  of  any  PCM  system, 
the  size  and  shape  of  the  PCM  container  must  correspond  to  the 
melting  time  of  the  PCM  and  the  daily  insolation  at  a  given 
location,  if  the  source  of  energy  is  a  solar  collector.  PCMs  are 
typically  placed  in  long  thin  heat  pipes  [102],  cylindrical  containers 
[100,103]  or  rectangular  containers  [66,91  ].  A  survey  of  previously 
published  papers  dealing  with  LHTES  reveals  that  two  geometries 
commonly  employed  as  PCM  containers  are  the  rectangular  and 
cylindrical  containers.  The  most  intensely  analysed  LHTES  unit  is 
the  shell  and  tube  system,  accounting  for  more  than  70%.  This  is 
probably  due  to  the  fact  that  most  engineering  systems  employ 
cylindrical  pipes  and  also  heat  loss  from  the  shell  and  tube  system 
is  minimal.  Fig.  1  gives  the  schematics  of  the  cylindrical  and 
rectangular  containers.  Table  4  lists  some  of  the  references  to  the 
studies  undertaken  employing  the  two  most  common  container 
geometries. 

3.2.  Configurations  of  cylindrical  PCM  containers 

Three  modes  of  cylindrical  PCM  container  configurations  are 
distinguished.  The  first  is  where  the  PCM  fills  the  shell  and  the  heat 
transfer  fluid  flows  through  a  single  tube  (Fig.  la)  [60,80,98-100] 
designated  the  pipe  model.  In  the  second  model  the  PCM  fills  the 
tube  and  the  HTF  flows  parallel  to  the  tube  (Fig.  lb)  [56,107,80]. 
According  to  Esen  et  al.  [80],  who  studied  the  two  models 
theoretically  by  comparing  the  effects  of  various  thermal  and 
geometric  parameters;  cylinder  radii,  total  PCM  volume,  mass  flow 
rates  and  inlet  temperatures  of  HTF  on  the  storage  time 
recommended  the  pipe  model  because  it  recorded  a  shorter  melt 
time.  This  was  because  the  thicker  the  PCM  mass,  the  longer  the 
melt  time  of  the  PCM.  An  added  advantage  for  the  pipe  model  not 
mentioned  is  the  fact  that  the  pipe  model  has  a  lower  heat  loss  rate 
to  the  environment  because  most  heat  supplied  from  the  center 
ends  up  heating  the  PCM.  The  third  cylinder  model  is  the  shell  and 
tube  system  [60,99]  commonly  used  to  improve  heat  transfer  in 
PCMs.  Agyenim  et  al.  [99]  conducted  an  experimental  energy 
storage  system  to  compare  horizontal  shell  and  tube  heat 
exchanger  (4  tubes)  and  a  pipe  model  incorporating  a  medium 


Table  2 

Thermophysical  properties  of  PCMs  investigated  for  different  applications. 


Compound 

Reference(s) 

Melting  temp,  Tm  (°C) 

Heat  of  fusion, 

Z  (kjkg-1) 

Specific  heat  capacity, 

Cp  (kj  kg"1  K"1) 

Thermal  conductivity, 
k  (Wm-1  K"1) 

Density,  p  (kgm  3) 

Water-ice 

[66-70] 

0 

335 

4.2 

2.4  (liquid)  0.6 

1000 

GR25 

[74] 

23.2-24.1 

45.3 

1.2  (solid)  1.2  (liquid) 

- 

- 

RT25-RT30 

[72,73] 

26.6 

232.0 

1.80  (liquid)  1.41  (solid) 

0.18  (liquid)  0.19  (solid) 

749  (liquid)  785  (solid) 

n-Octadecane 

[55,56,74-79] 

27.7 

243.5 

2.66  (liquid)  2.14 

0.148  (liquid)  0.190  (solid) 

785  (liquid)  865  (solid) 

CaCl2-6H20 

[66,80,81] 

29.9 

187 

2.2  (liquid)  1.4  (solid) 

0.53  (liquid)  1.09  (solid) 

1530  (liquid)  1710  (solid) 

Rubitherm  RT 

[37] 

Na2S0410H20 

[60,80,81] 

32,  39 

180 

2.0  (liquid)  2.0  (solid) 

0.15  (liquid)  0.3(solid) 

1460  (solid) 

Paraffin  wax 

[72,80,81] 

32-32.1 

251 

1.92(solid)  3.26(liquid) 

0.514  (solid)  0.224  (liquid) 

830 

Capric  acid 

[82,83] 

32 

152.7 

- 

0.153  (liquid) 

878  (liquid)  1004  (solid) 

Polyethelene  glycol  900  (PEG900) 

[84] 

34 

150.5 

2.26  (liquid)  2.26  (solid) 

0.188  (liquid)  0.188  (solid) 

1100  (liquid)  1200  (solid) 

Lauric-palmitic  acid  (69:31)  eutectic 

[85] 

35.2 

166.3 

2.41  (liquid)  1.77  (solid) 

Laurie  acid 

[82,85-87] 

41-43 

211.6 

2.27(liquid)  1.76(solid) 

1.6 

1.76(solid)  0.862  (liquid) 

Stearic  acid 

[61,87-89] 

41-43  (67-69)* 

211.6 

2.27  (liquid)  1.76  (solid) 

1.60  (solid) 

862  (liquid)  1007  (solid) 

Medicinal  paraffin 

[60] 

40-44 

146 

2.3  (liquid)  2.2  (solid) 

2.1  (liquid)  0.5  (solid) 

830  (solid) 

Paraffin  wax 

[90] 

40-53 

PI  16-Wax 

[55,60,80,81] 

46.7-50 

209 

2.89  (liquid)  2.89  (solid) 

0.277  (liquid)  0.140  (solid) 

786  (solid) 

Merck  P56-58 

[53] 

48.86-58.06 

250 

2.37  (liquid)  1.84  (solid) 

- 

- 

Commercial  paraffin  wax 

[94] 

52.1 

243.5 

- 

0.15 

809.5  (solid)  771  (liquid) 

Myristic  acid 

[35,58,95] 

52.2 

182.6 

- 

- 

- 

Paraffin  RT60/RT58 

[93,94] 

55  to  60 

214.4-232 

0.9 

0.2 

775  (liquid)  850  (solid) 

Palmitic  acid 

[83,85,87,89,92,95] 

57.8-61.8 

185.4 

0.162  (liquid) 

850  (liquid)  989  (solid) 

Mg(N03)2-6H20 

[71.88] 

89 

162.8 

- 

0.490  (liquid)  0.611  (solid) 

1550  (liquid)  1636  (solid) 

RT100 

[96] 

99 

168 

2.4  (liquid)  1.8  (solid) 

0.2  (liquid)  0.2  (solid) 

770  (liquid)  940  (solid) 

MgCl2-6H20 

[71,97] 

116.7 

168.6 

2.61  (liquid)  2.25  (solid) 

0.570  (liquid)  0.704  (solid) 

1450  (liquid)  1570  (solid) 

Erythritol 

[50,98-100] 

117.7 

339.8 

2.61(liquid)  2.25(solid) 

0.326  (liquid)  0.733  (solid) 

1300  (liquid)  1480  (solid) 

Na/K/NOs  (0.5/0.5) 

[54,65] 

220 

100.7 

1.35 

0.56 

1920 

ZnCl2/I<Cl  (0.319/0.681) 

[54,65] 

235 

198 

- 

0.8 

2480 

NaN03 

[54,65] 

310 

172 

1.82 

0.5 

2260 

kno3 

[54,65] 

330 

266 

1.22 

0.5 

2110 

NaOH 

[54,65] 

318 

165 

2.08 

0.92 

2100 

KOH 

[54,65] 

380 

149.7 

1.47 

0.5 

2044 

ZnCl2 

[54,65] 

280 

75 

0.74 

0.5 

2907 

LiF-CaF2  (80.5:19.5)  mixture 

[101] 

767 

816 

1770  (liquid)  1.770  (liquid) 

1.70  (liquid)  3.8  (solid) 

2390  (liquid)  2390  (solid) 

-  Not  available. 
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Table  3 

Target  application  areas  for  some  PCMs  studied. 


619 


Temperature 
range  (°C) 

0-65 °C 


80-120°C 

>150 


PCMs  studied/melting  temperature  (°C) 


Target  application  area  (rationale  behind  selection  of  PCM) 


Paraffins  (-3  to  64),  water/ice/0,  stearic  acid/41-43, 
n-octadecane/27.7 


Erythritol/1 17.7;  RT100  (99);  MgCl2-6H20  (116.7) 
NaN03/310,  KN03/330  NaOH/318,  KOH/380  ZnCl2/280 


Storage  for  domestic  heating/cooling  (see  list  in  Table  1).  Passive 
storage  in  bio-climatic  building/architecture.  Thermal  storage  of  solar 
energy.  Application  in  off-peak  electricity  for  cooling  and  heating. 
Protection  of  electrical  devices. 

Storage  for  the  hot-side  of  LiBr/H20  absorption  cooling  system  with 
generator  temperature  requirements  of  less  that  120°C  [100]. 

Storage  for  solar  power  plants  based  on  parabolic  trough  collectors 
and  direct  steam  generation. 


Ic~  Ip 


HTF 


(a)  Pipe  model  (b)  Cylinder  model  (c)  Shell  and  tube  model  (d)  Rectangular  or  slab  container 
Fig.  1.  Classification  of  commonly  used  PCM  containers  in  terms  of  the  geometry  and  configuration. 


temperature  phase  change  material  (erythritol)  with  a  melting 
point  of  117.7  °C.  The  thermal  characteristics  in  the  systems  were 
analysed  using  isothermal  contour  plots  and  temperature  time 
curves.  Temperature  gradients  along  the  three  directions  of  the 
two  systems;  axial,  radial  and  angular  directions  were  also 
analysed  and  compared.  Heat  transfer  in  the  shell  and  tube 
system  was  found  to  be  dominated  by  the  effect  of  multiple 
convective  heat  transfer  compared  to  conductive  heat  transfer  in 
the  pipe  model.  The  temperature  gradients  recorded  in  the  axial 
direction  for  both  the  pipe  and  shell  and  tube  systems  during  the 
change  of  phase  were  reported  to  be  2.5%  and  3.5%  that  of  the  radial 
direction  respectively,  indicating  essentially  a  two-dimensional 
heat  transfer  in  both  systems.  The  onset  of  natural  convection 
through  the  formation  of  multiple  convective  cells  in  the  shell  and 
tube  system  significantly  altered  the  shape  of  the  solid  liquid 
interface  fluid  flow  and  indicated  complete  melt  time  within  5  h 
compared  to  more  than  8  h  for  the  pipe  model.  The  authors 
recommended  the  shell  and  tube  system  for  the  charging  of  PCMs. 


Table  4 

References  to  container  geometry  employed  in  the  study  of  PCMs. 


Geometry 

Reference(s) 

Cylindrical  (concentric  annulus) 

[60,66,74,76,78,80,81,86,93,95, 

97,102,106,107] 

Cylindrical  (shell  and  tube) 

[56,75,90,101,104,105] 

Rectangular/slab 

[66,68,72,74,77,79,91,108] 

Spherical 

[74] 

3.2.  Counter-current  and  parallel  HTF  flow  directions 

In  a  cylindrical  container  assembly,  two  possibilities  exist  for 
the  flow  direction  of  the  heat  transfer  fluid  during  charging  and 
discharging  of  the  PCM  energy.  The  two  modes  are  either  the  hot 
and  cold  fluids  are  introduced  into  the  tube  from  the  same  end 
(parallel  flow)  or  the  hot  and  cold  fluids  are  introduced  from  the 
opposite  ends  (counter-current  flow)  of  the  heat  transfer  tube 
during  charging  and  discharging  respectively.  Fig.  2  illustrates  the 
schematic  diagram  of  the  parallel  and  counter-current  flow 
principles.  For  each  pair,  the  upper  arrow  represents  the  direction 
of  HTF  flow  during  charging  and  the  lower  arrow  represents 
discharge  direction  of  the  HTF. 

Results  from  the  numerical  simulations  conducted  by  Gong  and 
Mujumdar  [101]  to  investigate  the  effect  of  the  parallel  and 
counter-current  flow  modes  using  a  mixture  of  80.5%  LiF  and  19.5% 
CaF2  as  a  PCM  and  He/Xe  mixture  as  a  HTF  showed  that  parallel 
flow  increases  the  energy  charge/discharge  rate  by  5%  more  than 
counter-current  flow.  This  was  because  the  penetration  depth  of 
the  solid-liquid  phase  change  interface  during  the  charge/ 
discharge  was  larger  due  to  higher  temperature  difference  at 
the  fluid  inlet  if  hot  and  cold  fluids  enter  from  the  same  end.  In 
addition,  supercooling  of  the  PCM  did  not  occur  in  the  fluid  inlet 
region  and  heat  transfer  between  the  heat  transfer  fluid  and  the 
PCM  did  not  deteriorate.  On  the  other  hand,  counter-current  flow 
for  the  charge/discharge  processes  produced  significant  super¬ 
cooling  of  the  PCM  in  the  inlet  region  of  the  cold  fluid. 
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Charging 


Discharging 


(a)  Parallel  HTF  flow 


l 


Fig.  2.  The  physical  model  illustrating  parallel  and  counter-current  HTF  flow  in  a 
shell  and  tube  system. 


Charging 


Discharging 


3.3.  Investigating  parameters  in  PCM  containers 

Apart  from  the  container  geometry  and  configuration,  various 
thermal  and  geometric  parameters  are  known  to  affect  the  thermal 
performance  of  LHTES.  Table  5  lists  the  references  to  some  of  the 
operating  parameters  investigated  and  their  key  findings. 

Key  findings  from  numerical  predictions,  tested  against 
experimental  data  for  various  operating  parameters  have  shown 
that,  in  order  to  optimize  the  performance  of  a  PCM  storage  unit, 
the  thermal  and  geometric  parameters  listed  in  Table  5  must  be 
carefully  selected.  Literature  reviewed  has  shown  a  superior 
performance  using  the  shell  and  tube  configuration  followed  by 
the  pipe  model  with  the  PCM  at  the  shell  side  and  the  heat  transfer 
fluid  flowing  through  the  center.  It  has  therefore  been  recom¬ 
mended  by  many  authors  including  [60,75,99,118]. 

4.  Heat  transfer  in  PCMs  and  enhancement  techniques 

Most  PCMs  have  unacceptably  low  thermal  conductivity, 
leading  to  slow  charging  and  discharging  rates,  hence  heat  transfer 
enhancement  techniques  are  required  for  most  LHTES  applica¬ 
tions.  Several  studies  have  been  conducted  to  study  heat  transfer 
enhancement  techniques  in  phase  change  materials  (PCMs)  and 
include  finned  tubes  of  different  configurations  [59,69,94,97,98, 
102,109-112],  bubble  agitation  [112],  insertion  of  a  metal  matrix 
into  the  PCM  [113,114],  using  PCM  dispersed  with  high 
conductivity  particles  [115],  micro-encapsulation  of  the  PCM 
[116,117]  or  shell  and  tube  (multitubes)  [99,118].  Fig.  3  shows  the 
pictures  and  drawings  of  some  of  the  common  heat  transfer 
enhancement  techniques  studied.  Table  6  lists  some  examples 
where  enhancement  techniques  have  been  investigated. 

Majority  of  the  heat  enhancement  techniques  have  been  based 
on  the  application  of  fins  embedded  in  the  phase  change  material. 
This  is  probably  due  to  the  simplicity,  ease  in  fabrication  and  low 
cost  of  construction.  This  is  followed  by  the  impregnation  of  metal 


matrix  into  the  PCM  using  high  conductivity  materials  such  as 
carbon  fiber  and  brushes  and  multitubes.  A  selection  of  four  studies 
where  fins  (circular  and  longitudinal),  metal  matrix  and  shell  and 
tube  systems  approach  have  been  adapted  are  described  in  Table 
7.  It  can  be  seen  from  Table  7  that  different  approaches  were  used 
in  the  analysis  with  different  results  and  interpretations. 

The  general  observations  drawn  from  the  various  studies 
demonstrate  that,  irrespective  of  the  PCM  used,  the  heat  transfer 
characteristics  of  the  PCMs  can  be  improved  using  all  of  different 
enhancement  techniques.  Choi  and  Kim  [97],  Horbaniuc  et  al. 
[102],  Velraj  et  al.  [93]  and  Hamada  et  al.  [75],  employed 
different  experimental  setups,  different  container  materials,  and 
different  PCMs  to  investigate  the  heat  transfer  enhancement 
characteristics  of  PCMs. 

In  terms  of  performances  of  heat  transfer  enhancement 
techniques  and  systems  used  by  Choi  and  Kim  [97],  Horbaniuc 
et  al.  [102],  Velraj  et  al.  [93]  and  Hamada  et  al.  [75],  the  best 
enhancement  technique  as  reported  in  the  literature  was  that  due 
to  Velraj  et  al.  [93]  where  the  effective  thermal  conductivity 
calculated  employing  paraffin  with  lessing  rings  was  ten  times 
(2  Wrrr1  K_1)  greater  than  the  thermal  conductivity  of  paraffin 
(0.2  W  itt1  K_1).  This  must  be  qualified  by  the  fact  that  different 
researchers  employed  different  experimental  setups,  different 
container  materials,  and  different  PCMs  to  investigate  the  heat 
transfer  enhancement  characteristics  (Table  8).  For  example 
during  the  charging  and  discharging  stages  of  the  experiment,  the 
heat  transfer  fluid  in  the  experiment  from  Horbaniuc  et  al.  [102] 
passed  over  just  a  fraction  of  the  whole  setup  rendering  heat 
addition  and  removal  inefficient.  Results  derived  from  any  such 
setup  cannot  be  considered  to  be  optimal. 

Different  researchers  used  different  parameters  to  assess  the 
heat  transfer  enhancement  in  the  PCMs.  Velraj  et  al.  [93]  evaluated 
the  enhancement  of  the  heat  transfer  using  the  effective  thermal 
conductivity  taken  from  a  two-dimensional  enthalpy-temperature 
governing  equation  which  assumed  no  variation  of  temperature 
and  thermal  conductivity  in  the  axial  direction.  Results  were 
presented  graphically  using  temperature-time  curve  and  as  such 
limits  the  application  tending  to  other  applications.  Horbaniuc 
et  al.  [102]  measured  performances  of  fins  in  terms  of  the  interface 
freezing  stage  and  the  time  taken  for  complete  solidification  to  be 
achieved  using  parabolic  and  exponential  approximations. 
Hamada  et  al.  [75]  used  the  effective  thermal  conductivity 
proposed  by  Fukai  et  al.  [126]  to  assess  and  compare  results 
with  the  control  system  with  no  heat  transfer  enhancement.  In  the 
case  of  Choi  and  Kim  [97],  the  key  parameters  used  to  assess  the 
heat  transfer  enhancement  of  the  circular  finned  system  were  the 
ratio  of  overall  heat  transfer  coefficient  in  the  finned  and  the 
unfinned  tube  systems.  They  reported  a  ratio  of  3.5  for  a  surface 
area  ratio  of  3.2  between  the  finned  and  the  unfinned  tube  systems. 
A  comparison  was  also  made  by  deriving  a  relationship  for  the  ratio 
of  the  total  amounts  of  heat  recovered  in  the  finned  and  unfinned 
tube  systems  correlated  with  dimensionless  parameters  of  Fourier, 
Stefan  and  Reynolds  with  a  correlation  coefficient  of  0.994  and  a 


Table  5 

Examples  of  common  operating  parameters  investigated  in  the  study  of  phase  change  containers. 


Reference(s)  Parameters )  investigated 


[561 

1.  PCM  cylinder  (shell)  radii 

2.  Mass  flow  rate 

3.  HTF  inlet  temperature 

[661 

Comparison  of  rectangular  and 
cylindrical  containers 

[74] 

Melt  time 

[101] 

Effect  of  the  parallel  and 
counter-current  flow  modes 

Findings 

1 .  Increasing  mass  flow  rate  increases  the  heat  transfer  rate  and  the  molten  volume  fraction. 

2.  Increasing  inlet  temperature  rate  increases  the  heat  transfer  rate  and  the  molten  volume  fraction. 

For  larger  quantities  of  PCM  (>2kg),  the  melting  time  of  cylindrical  container  was  nearly  twice 
that  of  the  rectangular  one  for  equal  volume  and  equal  heat  transfer. 

Average  and  total  heat  fluxes  (4852  kj/m2  and  -4403  kj/m2)  approach  the  latent  heat  storage 
capacities  per  unit  surface  area  as  change  in  temperature  approaches  zero. 

Energy  charge/discharge  rate  for  parallel  flow  is  5%  more  than  counter-current  flow  mode 
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(i)  Longitudinal  or 
axial  fins 


(x)  Finned  Re  ctangular 
Container 


(xv)  Polypropylene 
flat  panel 


(ii)  Circular  fins 


(xi)  Graphite  Hakes 


(xvi)  Module  beam 


(iii)  Multitubes  or 
shell  and  tube 


(vii)  Encapsulation 


(xii)  Steel  metal 
ball  capsules 


(xvii)  PCM- 
Graphite 


(ix)  Metal  Matrix 


(xiv)  Polyolefine 
spherical  balls 


(xviii)  Compact  flat 
panel 


Fig.  3.  Heat  transfer  enhancement  methods  employed  in  phase  change  material  research. 


Table  6 

Heat  transfer  enhancement  techniques  employed  in  PCMs. 


Reference(s)  Nature  of  study 


Description  of  enhancement  technique 


[119] 

Experimental 

[120] 

Analytical 

[118] 

Numerical  and  experimental 

[121] 

Numerical 

[112] 

Numerical  and  experimental 

[122] 

Experimental 

[110] 

Numerical 

[79] 

Numerical  and  experimental 

[123] 

Analytical 

[124] 

Experimental 

[125] 

Numerical 

[99,100] 

Experimental 

LiH  PCM  encapsulated  in  smaller  container  and  placed  in  a  bigger  cylindrical  container  5%  Ni  added  uniformly  to  LiH. 
Aluminium  and  steel  internal  fins  in  the  solidification  of  n-octadecane  and  salt-hydrate  climsel. 

Staggered  cylindrical  tubes  (multitubes)  in  a  shell  and  tube  system  using  Indonesian  traditional  substance 
called  Mikro  as  a  PCM. 

Plate  type  fins  mounted  on  the  vertical  cooling  surface  of  a  rectangular  vessel  to  study  solidification  of  nitric 
molten  salt. 

Internal  longitudinal  fins  to  improve  heat  transfer  in  paraffin  PCM  in  a  cylindrical  container. 

A  spiral  twisted  copper  fins  in  a  cylindrical  tube  to  enhance  heat  transfer  during  the  melting  of  stearic  acid. 

Melting  of  ice  in  a  porous  medium  surrounded  by  finned  surfaces. 

System  of  7  aluminium  rectangular  containers  filled  with  n-octadecane  and  used  as  fins  to  improve 
heat  transfer  of  PCM. 

Simplified  analytical  model  to  study  heat  transfer  of  fins  and  honeycombs  in  PCMs. 

Ultrasonic  energy  used  to  enhance  heat  transfer  in  Erythritol  PCM. 

Micro-encapsulation  of  phase  change  slurries  to  improve  convective  heat  transfer  in  the  PCM. 

Internal  longitudinal  and  circular  fins  and  shell  and  tube  system  to  improve  heat  transfer  in  Erythritol. 
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Table  7 

A  comparison  of  results  of  study  using  fins,  metal  matrix  and  shell  and  tube  system. 


Reference(s) 


[97] 


[102] 


[93] 


[75] 


Heat  enhancement 
technique 

Circular  fins 


Longitudinal  fins 


Fins,  rings  and 
bubble  agitation 


Carbon  fiber 
and  brushes 


Method 


Analysis  and  results 


Freezing  of  0.36  kg  weight 
MgCl2-6H20 


Unfinned  tube  system:,  7=32—  =  2.39x 

J  Cimax 

10' -2Fo0J65Ste°- 528 Re0 195 ,  Re  <  7000 


Finned  tube  system:  q3L_ 
10 -3Fo°-781Ste°-616Re0-^ 


=  2.23x 
Re  > 7000 


Solidification  of  PCM 
(Tm  =  60°C) 


Parabolic  approximation  of  temperature 
profile: 


ATfin~AT2finad\n  2^/ATfin-Z2a(t) 
2  y/*Tfin-2Zy/m 


2 


Solidification  of  Paraffin 
RT  60  and  RT  58 


Melting  of  PCM 
n-octadecane 


Exponential  approximation  of  temperature 


profile:  =  akxt  1//4,  0 


AT(x,Q 


At 


fin 


and 


_  2  Xpsks 

k  k2finSfinATfin 

Recommended  internal  longitudinal  fins  and 
lessing  rings  for  solidification  enhancement 
and  the  use  of  an  evacuated  tube  for  melting 
when  using  paraffin. 


Estimation  of  effective  thermal  conductivities 
using  the  semi-empirical  equation  from  Fukai 


et  al. 


s ’ill 


—  0.456{1  -  (1  -Xflber)2/3} 


(k  fiber /kp)  +  (1  ^  fiber) 


1/3 


Observation  or  recommendation 


Radial  temperature  gradient  in 
the  finned  tube  system  is  uniformly 
larger  than  the  unfinned  tube  system. 
The  axial  temperature  difference  for 
the  finned  tube  was  smaller  than  for 
the  unfinned  tube  in  the  liquid  phase. 
Ten  fins  would  be  sufficient  to  allow 
a  complete  solidification  of  PCM  with 
Tm  =  60  °C  in  an  imposed  freezing 
time  of  180  min. 


Effective  thermal  conductivity  keff 
calculated  by  employing  paraffin 
with  lessing  rings  was  2Wm"loC_1 
(ten  times  greater  than  the  thermal 
conductivity  of  paraffin). 

Heat  brush/PCM  composite  was  found 
to  be  inferior  to  the  fiber-chip/PCM 
composite  in  terms  of  heat  transfer 
rate  of  the  bulk  but  the  former  was 
found  superior  to  the  latter  in  terms 
of  overall  heat  transfer  rate. 


Table  8 

Container  materials  used  by  previous  researchers  for  PCM  heat  transfer  enhancement. 


Reference(s) 

PCM  container  material  and  position 

Heat  transfer  tube 

Enhancement  technique 

[971 

Pyrex  glass/vertical 

Stainless  steel  tube 

5  stainless  steel  circular  fins 

[102] 

Horizontal 

- 

Longitudinal  fins 

[93] 

Aluminium/vertical 

Aluminium  tube 

Fins,  rings,  bubble  agitation 

[75] 

Acrylic  resin/vertical 

4  steel  tubes 

Carbon  fiber  and  brushes 

standard  deviation  of  0.023  and  0.028,  respectively.  The  above 
discussion  illustrate  the  fact  that  there  is  no  unified  international 
or  national  standard  methods  (such  as  British  Standards  or  EU 
standards)  developed  to  test  PCMs,  making  it  difficult  for 
comparison  to  be  made  to  assess  the  suitability  of  PCMs  to 
particular  applications.  A  standard  platform  such  as  British 
Standards,  EU  standards  needs  to  be  developed  to  ensure  same 
or  similar  procedure  and  analysis  (performance  curves)  to  allow 
comparison  and  knowledge  gained  from  one  test  to  be  applied  to 
another.  Also  contradictions  exist  in  the  thermophysical  properties 
of  PCMs  provided  especially  the  latent  heat  of  fusion,  thermal 
conductivity  and  densities  in  solid  and  liquid  states.  This  is  again 
due  to  the  absence  of  unified  certification  standards  and 
procedures. 

5.  Characterising  the  effect  of  melting  and  solidification  rates 
in  PCMs 

To  design  PCM  systems  requires  a  good  understanding  of  the 
fundamental  heat  transfer  processes  involved  in  accurately 
predicting  the  thermal  performance  of  the  PCM  system.  Based 
on  both  experimental  and  numerical  investigations  conducted  to 
examine  the  thermal  characteristics  of  LHTES  systems,  various 
correlations  have  been  developed  relating  thermal  performance 
and  dimensionless  numbers  in  given  parametric  domains. 
Dimensionless  numbers  are  widely  applied  in  modelling  the 
PCM  systems  in  order  that  the  knowledge  gained  from  one  study 
can  be  extended  beyond  its  source  of  acquisition.  Among  the 
thermal  performance  parameters  for  which  correlations  have  been 


developed  are  melted  volume  fraction,  temperature  profile,  melt 
time  and  melting  rate.  Table  9  gives  the  mathematical  definitions 
of  the  common  dimensionless  numbers  used  in  the  study  of  LHTES 
systems  and  their  significances  to  the  phase  change  process.  Table 
10  gives  references  to  some  of  the  correlations  derived  using 
dimensionless  numbers. 

In  spite  of  the  significant  literature  available  in  quantifying 
melting  and  solidification  rates,  there  has  been  no  attempt  to 
correlate  available  data.  The  two  main  reasons  for  this  are: 

•  individual  authors  used  different  phase  change  materials  with 
different  heat  transfer  characteristics.  In  the  case  where  the  same 
PCM  has  been  used,  the  researchers  employed  different 
dimensionless  parameter  ranges  and  presentations  making  it 
difficult  to  cross-correlate  between  the  characteristics  influen¬ 
cing  the  heat  transfer  in  specific  PCMs. 

•  some  researchers  presented  results  without  using  dimensionless 
parameters  making  it  impossible  to  extend  the  knowledge  to 
applications  beyond  the  original  source. 

For  example  the  empirical  correlation  for  the  melted  volume 
fraction  from  Wang  et  al.  [84],  for  example  predicted  the 
experimental  data  with  an  average  error  of  9.3%  and  the  deviation 
was  attributed  to  the  sensible  heating  that  takes  place  at  the  very 
early  times  due  to  the  presence  of  subcooling  within  the  PCM.  In 
the  case  of  the  time  averaged  Nusselt  number,  the  correlation  was 
found  to  fit  the  individual  experimental  data  with  an  average  error 
of  8.8%.  The  dependency  of  the  melted  volume  fraction  on  Fo  and 
Ste  but  minimal  for  Ra  was  that,  the  dimensionless  groups  Fo  and 
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Table  9 

Definitions  of  some  of  the  common  dimensionless  parameters  employed  in  the  analysis  of  phase  change  problems. 
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Number 

Source 

Formula 

Biot,  Bi 

[126] 

Bi  —  - 

dl  —  k 

Nusselt,  Nu 

[75,84] 

Nu  =  f 

Stefan,  Ste 

Dimensionless  time  or  Fourier  number,  Fo 

[127] 

[75] 

Ste  =  CpfT 

T  =  kt2 
pc  pi 

Rayleigh,  Ra 

[84] 

Ra  =  sfiAT'3 

av 

Prandtl,  Pr 

[84] 

Pr  =  % 

Reynolds  number,  Re 

[122] 

Re  =  V 

Grashof  number,  Gr 

[128] 

Cr  =  +"3 

VZ 

Significance  (determination) 

Ratio  of  conductive  to  convective  heat  transfers  resistance. 

Determines  uniformity  of  temperature  in  solid. 

Ratio  of  the  conductive  thermal  resistance  to  convective  thermal 
resistance.  Determines  the  ratio  of  actual  heat  transferred 
by  a  moving  fluid  to  the  heat  transfer  that  would  occur  by  conduction. 
Ratio  of  thermal  capacity  of  the  melted  solid  to  the  latent  heat. 
Characterises  heat  flux  into  a  body  or  system. 

Determines  the  onset  of  convection.  Below  a  critical  value, 
heat  transfer  is  primarily  conduction. 

Approximates  the  ratio  of  momentum  diffusivity  to  thermal 
diffusivity.  Low  Pr  means  effective  heat  conduction  with 
dominant  thermal  diffusivity.  High  Pr  means  effective  heat 
convection  with  dominant  momentum  diffusivity. 

Ratio  of  inertial  forces  to  viscous  forces.  Determines 
whether  flow  is  laminar  or  trubulent. 

Approximates  the  ratio  of  buoyancy  force  to  the  viscous  force. 


Table  10 

A  list  of  some  of  the  correlations  derived  by  different  references  using  dimensionless  numbers. 


Reference(s) 


Correlation  derived 


[84] 


[74] 


[78] 


v~0  =  a(t-tc ) 

X-  =  4.73Fo°-906Ste1-538Ra0002 

Vn 


JVii  =  0.219Ra°-387  Pr0019(f) 


0.062 


tf= _ l- _ 

Lf  2a  (1  +w)Ste 


{1  +(0.25  + 0.1 7w°-7)5te} 


q  =  =^{1  +  (0.1210  +  0.0424w)Ste(°-7645”°-2022w) } 


keq(  t)  =  0.2280  Ra(r) 


1/4 


1  _  d(r)" 

1  d 


1/4 


[108] 

[97] 


A(t)[1  -  lnA(r)]  =  1  -  4  J0T  keq(T)dr 


M(f)  =  l/o1  /20(y,  t)dy,  He(f)  =  /0‘  \dy  C(f)  =  /0‘  fdy  with  M  =  He-C 
^  =  2.39  x  10 -2Fo0J65Ste0-528  Re0195 

Li  max 


-1  2  A,-, 


1  2  B  rhl 


=  2.23  X  lCr3Fo0-781Ste°-616Re0-503 

Limax 


Interpretation  of  correlation 

Linear  function  of  the  melted  volume  of  PCM  with  time. 

Melted  volume  fraction  for  the  entire  duration  of  test  runs  in 
terms  of  Fo,  Ste  and  Ra. 

Time  averaged  Nusselt  number  for  the  entire  melting  process. 

Total  heat  input  per  unit  surface  area  Q  and  average  heat  flux 
at  r=l  over  the  duration  of  the  melting  process  tf. 

Equivalent  thermal  conductivity  and  dimensionless 
melting  rate  [1  -A(t)]. 


Instantaneous  melting  rate,  instantaneous  heating  rate  and 
instantaneous  cooling. 

Radial  temperature  gradient  in  the  unfinned  and  finned  tube 
PCM  (MgCl2-6H20)  system. 


Ste  were  the  important  factors  to  consider  in  the  phase  change 
process. 

The  difficulty  in  comparing  data  was  illustrated  by  Agyenim 
[96]  by  comparing  the  dependence  of  Nusselt  number  on  Reynolds 
number  for  four  independent  investigators: 


i.  Hamada  et  al.  [75]  employed  equation  from  Petukhov  [129], 
Gnielinski  [130]. 


Nu  = 


( / /nc/2)  (Re  -  1 000)Pr 


l  +  12.7^//ric/2(Pr2/3-l) 

(3000  <  Re  <  1 0s ,  0.5  <  Pr  <  2000)(  1 ) 


ii.  Zhang  and  Faghri  [111]  employed  empirical  correlation 
proposed  by  Edwards  and  Jensen  [131]. 


Nu 


0.023  (1/di)02 
(1.35  -0.351/di) 


ReosPr0.4 


with  l/di  calculated  from  available  data  as  0.50408. 

iii.  Sexana  et  al.  [132]  calculated  the  Nusselt  number  using 
equation  3.34  for  turbulent  flow. 

Nu  =  0.026 Pe08Pr1/3  (3) 

iv.  Esen  et  al.  [80]  and  Rohsenow  et  al.  [133]  employed  Dittus- 
Boelter  equation  for  turbulent  flow: 


Nu  =  0.023  Re0  8Pr°  4  for  (Re  >  2200) 


Agyenim  [96]  illustrated  that  although  all  correlations  showed  the 
same  general  trend  with  an  apparent  slight  differences  in  values 
recorded  for  the  lower  Reynolds  numbers  in  the  turbulent  regime, 
major  deviations  increased  with  increased  Nusselt  and  Reynolds 
numbers.  Independent  correlations  presented  demonstrated  that  in 
the  turbulent  region,  the  Nusselt  number  would  be  governed  by  the 
dimensionless  parameters,  Nu  =  Nu(Re,  Pr)  but  data  correlated  only 
marginally  well. 

5.2.  Phase  change  problem  formulation 

The  analysis  of  heat  transfer  problems  in  phase  change 
processes  is  complex  due  to  the  fact  that  the  solid-liquid 
boundary  moves  depending  on  the  speed  at  which  the 
latent  heat  is  absorbed  or  lost  at  the  boundary.  The  position 
of  the  boundary  is  unknown  and  forms  part  of  the  solution  Li  et 
al.  [134].  An  increasing  number  of  studies  have  been  conducted 
and  literatures  on  LHTES  problems  are  numerous.  Solutions  to 
phase  change  problems  include  analytical,  experimental  and 
numerical  using  one-dimensional,  two-dimensional  or  three- 
dimensional  models  to  solve  energy  formulated  equation. 
Energy  equation  is  formulated  in  various  ways  with  the  phase 
change  being  accounted  for  in  different  representations. 
Examples  of  some  of  the  formulation  by  various  references 
are  given  in  Table  11.  Table  12  provides  a  list  of  references  for 
the  various  solution  methods  employed  in  solving  the  phase 
change  problems. 
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Table  11 

Some  examples  of  formulation  for  the  phase  change  problem. 


Reference(s)  Governing  equation 


[74] 

[56] 

[101] 

[108] 

[121,135] 

[112] 


S  =  wir  (rW  S)  for  nt)  <r  <l,T  =  Tcr,  r  <  r(t ) 

{0  for  a  PCM  slab  insulated  at  one  end 
1  for  a  PCM  cylinder 
2  for  a  PCM  sphere 


f  =  ti  (“ri )  + 1  (“  i)  -  PAhf  w 

pCpicrf  g  =  27 zTj[/(r  -  Th)  -  rhCp  f 

Phase  change:  H(T )  =  h(T)  +  psf  Ah  f  where  h(T)  =  pCpdT 


PhC„.h  (^  +  l’%f)=f(Tp~TH)  +  kH 


Solid  phase  of  PCM:  +  Uo  =  ccs  |^r,  liquid  phase  of  PCM 


ar  _ 


dll; 


Mass  conservation  equation:  =  0 

Momentum  conservation  equations:  pi^  +  Uj^j  =^f- 


Conservation  of  energy:  pc  (  g  +  u  ,•  #  )  =  &  {  kM-, 


dT 


d  T 


dt 


lj  dX; 


dXj 


dX; 


:a  +  |  =  Ot/||  +  rf(T--T/)  =  0,f  +  |(ur)  +  |()  =  «/g 
PSfKT  -  Tref)  t  where  <r„  =  -  ^  (|J  +  £f)  and  p'  +  pgx2 


n  M  —  1 JL  (kr<yr\  I  1  _9_  fk  dT\ 

P  dt  —  r  dr  V,u  3r/  r  36>  Vr  36>J 

Phase  change  formula:  solid  region  (T< Tm - e):  H  =  cpT,  interface  (Tm-e<T<Tm  +  e):  H  =  cpT  +  (^)T  -  Tm  +  e,  liquid  region  (T> Tm  +  e):  H  =  cpT+h 


Table  12 

Numerical,  experimental  and  analytical  solution  and  validation  of  phase  change  problems. 


Reference(s)  Solution  and  validation 


Nature 


[66] 

[66] 

[80] 

[74] 

[104] 
[97] 
[93] 
[102] 

[75] 

[105] 
[81] 
[72] 
[68] 
[86] 
[91] 
[95] 

[76] 

[106] 

[77] 

[107] 

[78] 

[108] 

[79] 
[101] 
[60] 
[90] 


2D  fixed-grid  conjugate  enthalpy  based  method.  Finite  difference  and  tridiagonal  matrix  algorithm  solver. 
ID  enthalpy  based  method.  Gauss-Seidel  iteration  procedure. 

2D  conjugate  heat  transfer  enthalpy  based  method.  Gauss-Seidel  iteration  procedure 
ID  enthalpy  based  method.  Solved  via  Crank-Nicholson  scheme. 

ID  enthalpy  based  method.  Finite  difference  program-CYCLOP. 

ID  and  2D  enthalpy  method.  Implicit  finite  difference  method. 

1 -Parabolic  and  exponential  approximations  for  temperature  profile.  Newton-Raphson  technique. 

2D  enthalpy  method.  Solved  using  the  control  volume  method  by  Pantakar  (1980). 

2D  fixed-grid  conjugate  enthalpy  based  method.  Control  finite  difference  method. 

2D  enthalpy  based  method.  Gauss-Seidel  iteration  procedure  and  TRNSYS. 

2D  Navier-Stokes  and  energy  equations.  Finite  element  model. 

ID  and  2D  fixed-grid  enthalpy  based  method.  Control  volume  approach  by  Pantankar  (1980). 

ID  enthalpy  based  method. 

Fixed-grid  enthalpy-porosity  model.  Streamline  upwind/Petrov  Galerkin  finite  element  method. 

2D  fixed-grid  enthalpy  method.  Finite  volume  approach. 

2D  enthalpy  method.  Solved  using  Ozisik. 

ID  enthalpy  method. 

2D  mass,  momentum  and  energy  equations.  Finite  difference  method. 

ID  and  2D  enthalpy  method.  Fully  implicit  finite  difference  method. 

2D  energy  based  equation.  Finite  element  analysis. 

2D  fixed-grid  enthalpy  based  method.  Gauss-Seidel  iteration  method. 

3D  enthalpy  based  method.  Control  volume  method  by  Patankar. 


Numerical  and  experimental. 
Numerical  and  experimental. 
Numerical. 

Numerical. 

Analytical  and  Numerical. 
Experimental. 

Numerical  and  experimental. 
Analytical. 

Numerical  and  experimental. 

Numerical 

Numerical 

Numerical  and  experimental. 
Numerical. 

Experimental. 

Experimental. 

Experimental. 

Numerical. 

Numerical 

Theoretical 

Analytical. 

Analytical  and  experimental. 
Theoretical  and  experimental. 
Numerical 
Numerical. 

Theoretical. 

Experimental  and  numerical. 


A  survey  of  phase  change  formulation  reveals  that  the  most 
common  approach  applied  in  the  solution  of  phase  change 
problems  is  the  formulation  using  the  enthalpy  method.  The 
advantage  of  employing  the  enthalpy  procedure  has  been  reported 
as  the  removal  of  the  need  to  satisfy  conditions  at  the  phase  change 
front  which  most  authors  indicated  facilitate  the  implementation 
of  the  numerical  algorithm.  The  survey  of  publications  from  LHTES 
systems  with  cylindrical  containers  also  shows  that  nearly  half 
(47%)  of  the  studies  undertaken  are  numerical  and  half  of  the 
references  employed  two-dimensional  models  (Fig.  4). 

Results  from  the  numerical  methods  reported  appear  to  show 
that  they  offer  a  good  approach  to  solving  the  phase  change  problem 
although  most  of  the  available  solutions  to  phase  change  problems 
apply  to  one-  or  two-dimensional  systems  due  to  the  complexity  of 
the  equations  involved  in  the  phase  change.  The  latent  heat  phase 


change  problem  has  been  assumed  to  be  conduction  controlled  by 
some  authors  whereas  natural  convection  in  the  liquid  phase  has 
been  included  in  the  analysis  by  other  authors,  adding  to  the 
complexity  of  the  phase  change  problem.  Few  examples  of  phase 
change  problem  classification  have  been  reported  in  Table  13. 

The  numerical  formulation  has  also  concentrated  on  resolving 
specific  phase  change  material  problems  and  to  study  the 
characteristics  of  proposed  or  new  materials.  They  include: 

•  heat  transfer  during  melting  and  solidification  of  PCM;  evalua¬ 
tion  of  the  temperature  distribution  inside  PCM,  interface 
transient  location,  and  the  effect  of  natural  convection  in  the 
molten  zone  and, 

•  the  overall  thermal  performance,  system  heat  capacity,  charging 
and  discharging  rates  and  storage  temperatures. 


F.  Agyenim  et  al  /  Renewable  and  Sustainable  Energy  Reviews  14  (2010)  615-628 


625 


Numerical 

47% 


n.a  -  not 
available 


Fig.  4.  Percentages  of  PCM  solutions/analysis  by  methods  and  dimensions  for  28 
LHTES  systems  with  cylindrical  containers. 


Table  13 

PCM  heat  transfer  control  classification  by  different  references. 


Reference(s) 

PCM 

Control  classification 

[68] 

Water-ice  system 

Conduction 

[95] 

95%  purity  palmitic  acid 

Conduction 

[86] 

95%  purity  lauric  acid 

Conduction 

[76[ 

n-Octadecane 

Convection 

[106] 

- 

Convection 

[77] 

n-Octadecane 

Conduction/convection 

[107] 

Paraffin  wax,  paraffin  (p-116) 
and  stearic  acid 

Conduction/convection 

[78] 

n-Octadecane 

Conduction/convection 

Table  14 

Summary  of  method  of  analysis  with  key  findings  for  conduction/convection  control  phase  change  problems. 


Reference(s) 

PCM/method 

Parameters  used  to  examine  results 

Key  findings 

[86] 

Lauric  acid.  Pipe  model  with  the  PCM  on  the 
shell  side  and  the  heat  transfer  fluid  flowing 
through  the  center  of  the  tube. 

Reynolds  and  Stefan  numbers.  Graphical 
presentation  of  temperature  variations  in 
radial  and  axial  directions  for  different 
inlet  HTF  temperatures. 

Charging  time  was  less  than  solidification 
time  due  to:  (i)  melting  rate  being  governed 
by  convective  heat  transfer  and  (ii) 
conduction  thermal  resistance  for  the 
solidification  process  (Rcd  =  rH  ln(r/rw)/2/<s). 

[78] 

n-Octadecane.  PCM  inside  an  isothermally 
heated  horizontal  cylinder. 

Ra  and  Pr  numbers.  Photographic 
samples  at  dimensionless  times,  t. 

Graphical  presentation  of  equivalent 
thermal  conductivity  equation  derived  to 
estimate  dimensionless  melting. 

Melting  proceeded  faster  than  that  of  the  pure 
heat  conduction  and  heat  transfer  at  the  upper 
region  of  the  PCM  was  larger  than  the  lower 
region  due  to  the  presence  of  natural 
convection.  However,  decreasing  Ra  decreased 
convection  and  was  negligible  at  Ra<  105. 

[108] 

n-Octadecane.  PCM  enclosed  in  a  rectangular 
container  and  heated  at  a  constant  rate  from 
one  side. 

Nu,  thermal  conductivity  ratio  A,  solid 
subcooling  B  and  Ste,  Heating  rate  He 
and  Melting  rate,  M. 

Graphical  presentation  of  Nu,  A,  B,  M  and  He. 

A  has  practically  no  effect  on  the  heating  rate 

He  during  the  convection  dominated  regime 
and  has  a  minor  impact  on  He  during  the 
conduction  dominated  regime.  Higher  B  values 
produce  lower  melting  velocities,  especially  in 
the  early  stage  of  the  process.  Melting  rate 
decreases  as  Ste  increases. 

5.2.  Classification  of  heat  transfer  characteristics 

Problems  that  arise  due  to  heat  transfer  accompanied  by  a 
change  of  phase  can  be  classified  as: 

a.  conduction  controlled  phase  change, 

b.  convection  controlled  phase  change  or 

c.  conduction/convection  controlled  phase  change. 

References  to  some  of  the  control  classifications  for  some  of  the 
various  studies  reviewed  are  given  in  Table  13. 

Studies  based  on  the  assumption  that  conduction  is  the  major 
mechanism  of  heat  transfer  through  the  PCM  assumed  negligible 
convection  heat  transfer  in  the  melt.  In  the  case  of  convection 
controlled  phase  change,  the  common  assumption  was  that  the 
PCM  was  at  the  melting  temperature.  Models  based  on  pure 
conduction  and  pure  convection  have  not  been  able  to  predict 
accurately  the  melting  rates  nor  track  the  complex  motion  of  the 
solid-liquid  interface  at  all  times  during  the  phase  change.  Less 
than  10%  of  articles  surveyed  treated  convection  and  conduction  at 
different  stages  and  concluded  that  it  depends  on  the  phase 
transition  stage  to  determine  whether  or  not  conduction  or 
convection  controls  the  heat  transfer.  Among  these  are  Sari  and 
Kaygusuz  [86],  Hirata  and  Nishida  [78],  Agyenim  et  al.  [99]  and 
Zhang  and  Bejan  [108]. 

The  authors  accounted  for  both  conduction  and  convection, 
making  specifications  at  which  times  conduction/convection  is 
(are)  dominant  during  the  phase  change  process.  The  details  of  the 
work  by  Zhang  and  Bejan  [108]  together  with  Sari  and  Kaygusuz 
[86],  Agyenim  et  al.  [99]  and  Hirata  and  Nishida  [78]  establish  the 
following  stages  for  the  melting  process: 


a.  Sensible  heat  provision  at  the  start-up  period  during  which  heat 
added  increases  the  temperature  of  the  solid  PCM.  Pure 
conduction  heat  transfer  occurs  at  this  stage. 

b.  A  second  conduction  regime,  in  which  heat  is  purely  transferred 
by  conduction  from  the  heated  wall  to  the  PCM  and  solid-liquid 
interface  when  melting  had  just  begun. 

c.  The  transition  from  conduction  to  natural  convection  which 
starts  as  the  thickness  of  the  melt  layer  increases  and  the 
interface  starts  to  incline.  At  the  interface,  there  exists  at 
equilibrium,  a  solid  body  and  a  pool  of  its  own  liquid. 

d.  The  convection  regime,  when  most  of  the  solid  have  been 
melted  and  the  liquid  core  temperature  distribution  depends  on 
height  and  depth  and  not  on  time  (caused  by  buoyancy). 

Phase  change  problems,  first  treated  as  pure  conduction  controlled, 
has  in  recent  times  moved  to  a  different  level  of  complexity  with 
added  convection  in  the  melt  being  accounted  for  (Table  14). 

6.  Conclusion 

Studies  of  phase  change  systems  have  over  the  past  three 
decades  investigated  design  fundamentals,  system  and  process 
optimization,  transient  behaviour,  and  field  performance.  The 
research  and  development  has  been  broad  based  and  productive, 
concentrating  on  both  the  resolution  of  specific  phase  change 
materials  and  problems  and  the  study  of  the  characteristics  of  new 
materials.  The  review  reveals  the  different  heat  transfer  solution 
methods  employed  by  different  researchers.  They  include  theore¬ 
tical,  experimental  and  numerical  studies  that  have  been 
conducted  on  the  thermophysical  properties  of  new  phase  change 
materials. 
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The  following  conclusions  can  be  drawn  from  the  models 

reported  in  the  heat  transfer  and  phase  change  problems  reviewed: 

•  Most  of  the  research  on  phase  change  problems  has  been  carried 
out  within  the  temperature  range  0-60  °C  suitable  for  domestic 
heating/cooling.  Common  PCMs  studied  are  paraffins,  Glauber’s 
salts,  wax,  stearic  acid  and  n-octadecane.  Studies  of  medium 
temperature  latent  heat  thermal  energy  systems  to  be  applied  for 
hot-side  storage  of  absorption  air  conditioning  are  few. 

•  Most  phase  change  materials  have  low  thermal  conductivities 
requiring  heat  transfer  enhancement  techniques  to  improve 
their  rates  of  charging  and  discharging  of  energy.  Majority  of  the 
studies  undertaken  on  phase  change  materials  have  reported 
poor  thermal  conductivities  requiring  heat  transfer  enhance¬ 
ment  techniques.  As  such  current  studies  have  concentrated  on 
methods  of  heat  transfer  enhancement  using  mostly  fins, 
insertion/dispersion  of  high  thermal  conductivity  materials, 
multitubes  and  micro-  or  macro-encapsulation 

•  Numerical  methods  appear  to  offer  a  good  approach  to  solving 
the  phase  change  problem  although  most  of  the  available 
solutions  to  phase  change  problems  apply  to  one-  or  two- 
dimensional  systems  due  to  the  complexity  of  the  equations 
involved  in  the  phase  change. 

•  The  latent  heat  phase  change  problem  was  first  assumed  to  be 
conduction  controlled.  In  recent  times  however,  natural  con¬ 
vection  in  the  liquid  phase  has  been  included  in  the  analysis  of 
the  phase  change  problem  adding  to  the  complexity  of  the  phase 
change  problem. 

•  A  common  approach  observed  in  the  solution  of  phase  change 
problems  has  been  the  use  of  the  enthalpy  formulation.  The 
advantage  of  the  enthalpy  procedure  has  been  reported  as  the 
removal  of  the  need  to  satisfy  conditions  at  the  phase  change 
front  to  facilitate  the  implementation  of  the  numerical  algo¬ 
rithm. 

•  Variations  are  found  to  exist  in  correlations  derived  to 
characterise  the  effect  of  melting  and  solidification  of  PCMs. 
The  variations  have  been  attributed  to  the  different  methods  of 
analysis  and  the  choice  of  characteristic  lengths  used  and  as  such 
calls  for  caution  in  the  selection  of  correlation  for  design 
purposes  and  quantitative  comparisons. 

•  There  are  no  national  and  international  standards  for  testing 
PCMs.  As  such  researchers  have  over  the  years  applied  different 
approaches  in  both  test  procedures  and  analysis  of  results.  The 
development  of  a  unified  international  standard  for  test  and 
analysis  will  be  helpful  in  the  development  of  latent  heat  thermal 
energy  storage  systems.  Also  contradictions  exist  in  the 
thermophysical  properties  of  PCMs  provided  especially  the 
latent  heat  of  fusion,  thermal  conductivity  and  densities  in  solid 
and  liquid  states.  This  is  again  due  to  the  absence  of  unified 
certification  standards  and  procedures. 
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